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Abstract 

The Hubbard model is used to study an electronic system. In this paper we present 
the new path integral representation for Hubbard model. We have constructed the new 
supercoherent state for spinless electrons which appears from a set of eigenfunctions of 
atomic limit of strongly correlated systems. Exact calculation of nonlinear representa- 
tion of a supergroup has been carried. This group defines the transformation of atomic 
base. The general formalism we elaborate for Hubbard model is the one widely used 
in the gauge field theory of the nonlinear representation of a superconformal group. 
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1 Introduction 



The Hubbard model was originally constructed to describe a metal-insulator transition 
for spin-dependent fermions in a simple way [H El [3] . 

Today this model is still remain the main workspace for investigation of the strong elec- 
tron's correlation. There exist many approaches to this model for describing many electrons 
system: the band limit approximation for the weak interaction between electrons and the 
atomic limit for electrons with the strong coulomb repulsion. We start the series of papers 
in which we are intending to elaborate the new approach to the Hubbard model. 

We will present a new path integral approarch to the strong interaction regime. The 
main ingredients for us is the usage of a supercoherent state acting upon an atomic base and 
the work with effective functional for an electronic system. We will develop the procedure 
of geometric quantization for strongly correlated systems of electrons. Let us make the brief 
sketch of the program firstly developed in the series of papers of one of the authors [1]. 
One of the main distinction of this approach to the Hubbard model is that we treat the 
local supergroup of the local space-time as the main object of our theory. This supergroup 
generates the transformation of the space-time coordinates which assigns the arguments for 
any function describing this system. This representation of a supergroup in the superspinor 
space has to contain a Lorents or an SO (4) generators equal to the even subalgebra of the 
Hubbard operators. Odd Hubbard operators produce some superextension of the Lorents 
group into some supergroup. This dynamic supergroup is given by local dynamic superfields 
describing the local degrees of freedom of the strongly correlated electrons system. We shall 
introduce the supercoherent state depending on generalized angles equal to the bosonic and 
fermionic fields of the system. Parametrised by x,y,z,t space-time manifold which deter- 
mines the arguments of wave function defines us the spinor and superspinor bundle. This 
superbundle is determined by the supercoherent state. 

There are two kinds of the gauge fields in superspinor bundle: one sort of fields is the com- 
posite fields equal to the quadratic combination of odd grassmann fields and other sort are 
nonlinear fields determining the local coordinates frame of the four-dimensional space-time. 
In general, the local superspinor bundle defines the nonlinear representation of a super- 
conformal group as a maximal group of 4-dimentional space-time of interacting fermionic 
system. 

As the first step of quantization of electronic system with strong coulomb repulsion 
we perform the reformulation of Hubbard model into the atomic limit formalism. This 
approach is well known in Mott-Hubbard insulators theory. We want to point out that our 
approach includes all elements of geometric quantisation [S]: for example possess the algebra 
of 4-dimensional rotations (Lorents) group, Cartan differential one-forms which give us the 
lagrangian of system, the nonlinear representation of underlying supergroup as a ground for 
the supercoherent state. 

In the second paper of this series we will fulfill the exact calculation of the nonlinear 
representation of spinless supergroup generators of which appears in the so called "tower 
of symmetry" in the Hubbard model [1]. This gives us the possibility to introduce the 
supercoherent state for some deformed version of spinless algebra of the strongly correlated 
electron system and as a result to obtain the effective action in a future paper. 

Let us give the brief description of the calculation method for the finding representation 
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of dynamic supergroup in strongly interacting models. Our task is to find in this model such 
group structure which could help us to describe the specificity of strong correlation. We take 
the following construction as a base: 

1) we will collect space coordinates together with the time coordinate and will consider 
some curved space-time as a base in which Lorents subalgebra of superconformal group act 
in the spinor base on dynamical fields. 

2) full bases of electronic operators gives us the supergroup which is parametrised by 2 
dynamical fermionic fields: this fields comprise the real massless spinor. 

3) the superspinor representation of a supergroup gives us the supercoherent state de- 
scribed by the nonlinear function over odd grassmanian fields. This function characterises 
the local properties of the strongly correlated system. 

2 Atomic description of Hubbard model 

We consider the Hubbard model: 

ijcr i,a a,i 

here a^rja^j -electron creation and annihilation operators, n^^i - electron density 
operator W, U, -band width, one- site electron repulsion and chemical potential. 

At first we represent the Hubbard model in atomic bases which determine the atomic 
limit. This limit appear as a result of the following procedure. A zero approximation of the 
atomic limit is described by one-site repulsion term: 

i,a (T,i 

This hamiltonian can be diagonalized by the following one-site atomic eigenf unction: 

|0 >-; 1+ >-= tt^|0 y; |- ^= a||0 ^; |2 ^= «+a+|0 >- . (2) 

This bases gives us the fundamental representation of some supergroup in the space of 
dimension (2,2). Point out that the states |-|- |— are fermions but |0 |2 >- 
are bosons in our construction. We insert later odd grassmann fields and make this fact 
obvious ( ie states |-|- |— >-; will be depend on odd order of the grassmann fields but 
states |0 \'2 >- -on even order of the grassmann fields). All operators in this bases will 
be the matrix which are determined by commutation and anticommutation relation giving 
some superalgebra. Full set of the Hubbard operators have 16 operators part of which 

x°-, x+°, x^°, x-\ x^-) 

are the fermionic operators, but other part 

X++ - X—, ^2°, - X^^) 

-are the bosonic operators. X^^ -Hubbard operators contain only one non-zero element equal 
1 sitting on site in the matrix representation. Point out that this set of operators gives 
some bases for some superalgebra. 
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We have the following representation for creation-annihilating operators in this bases: 



a. 



+ 



The Hubbard model in this representation has the form: 



(3) 



(4) 



3 Supercoherent state for Hubbard model 

In constracting the supercoherent state we use the following interesting observation in in- 
terpretation of the set of atomic operators and function for on-site Hubbard repulsion. This 
observation can be formulated as the following statement: six even Hubbard operators con- 
stitute the subalgebra isomorphic with algera of Lorentz group or algebra of four dimensional 
rotation group in spinor representation. Complete derivation of this statement will be obtain 
in subsequent paper. 

To characterise the state of the system by coherent state we input some fields which 
depend on coordinates z and the time t. We have three component dynamic vector of 
the electricl field 

E = {E+{x, y, z, t), E-{x, y, zt), E'{x, y, z, t)), 
three component dynamical vector of the magnetic field 

h = {h^{x, y, z, t), h'{x, y, z, t), /i^(a;, y, z, t)) 

and four component dynamical odd grassmann fields 

X*{x,y,z,t),x{x,y,z,t), 

which are the fermionic fields giving the components of maiorana spinor. All dynamical 
fields appear in supercoherent state in the foUowng manner: 



G >— exp 



0>, 



(5) 



X* h h+ X 
X* h- -h, X 
E- X* X* -E, 

Exponent here act in space of atomic eigenfunctions (|0>,|-|->,| — >,|2>), Function 
I > is highest weight vector of representation space of supergroup which is given by 
exponent. 



4 Evolution operator for electronic system 

The transition amplitude of the evolution operator of the quantum systems is given by 
the following expresion: < Z/|e~*^'^*J^~*^)|Zj >. We want to obtain the expression for the 
effective functional using the states \Z >. Time evolution of the system is given by the 
following operator: 
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U{t,to) =Tordexp{-i H{T)dT) 

J to 

if t — to — St is small, ie St « 1, then 



rto+St 

U{to + St,to)^l-i / H{r)dr 

J to 



It is follow from this expression that the symbol for evolutionary operator has the fol- 
lowing form: 

rto+5t 

U{Z, Z*\to + St, to) = exp(-i / H(Z, Z*\r)dT). 

J to 

We devide time interval [to,t] by the number N and obtain N small intervals for finding 
the expression for symbol U {Z, Z*\t, to)- Consider the matrix elements of evolution operator 
exp{—iH{tf—ti)) between the states < Zf \ and \Zi >. Factorising operator exp{—iH{tf—ti)) 
by inserting the identity operator J d^{Z)\Z >< Z\ — 1 we obtain the following representa- 
tion: 

N 

< Zf\exp{-iH{tf - ti))\Zi >= / < Zf\ZN > 

k=i 

< Z^|e-^^^|Z^_i > .... < Zk-i\e-''"\Zy, > ... < Z^\e-''"\Zi >, 

here e = -^^jv^- In first order of e we can transform this formular and place the symbol 
of operator H in the exponent 

< Z,,,\e---\Z, > ^ <Z,,,\il-^eH)\Z,> ^ ^-.^g^ + o{e% (6) 
< Zk^i\Zk > < Zk^i\Zk > 

As a result we obtain the representation 



< Z/|e-^^(*^-*^)|Z, >= lim / rr dfi{Zk) < Zk+i\Zk > e ^'==^ <-k+i\H> , (7) 

N—>oo I 

fe=l 

here \Zo >= \Zi >; < Z^+il =< Zf\. Let define a variation of the following type \Z >: 
\SZk^i >— > —\Zk > . We have: 

N 

< Z/|e-*^(*/-*^)|Z, lim / nM/^(^ifc) < Zk\Zk >] < Zf\ZN > 

N—>oo I 

k=l 



exp{y{Ln{l ) - le ). 

in linear-slice approximation. We take the following expresion for the time derivative: 
Considering the first order in e we can take the folUowing expresion for the time derivative: 
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d\Z > _ \6Z > 
dt e 

In first order in e, we obtain the final path integral representation of the evolutionary 
operator in coherent state formalizm 



r\Z{tf)>=\Zf> 
J\Z{ti)>=\Zi> 



S[Z, Z* 



[ ^ dt [ dr{- 
Jti Jv 



< Z{r. f)\i^ - H\Z{r. f > 



(8) 



< Z{r,t\Z{r,t> 
-i[Ln(< Zf\Z(tf) >) - Ln(< Zi\Z{U) >)]). 
Measure of integration is given by the following expresion: 

D[Z,Z*]^ W \{dix[Z{r,t)\Z{T,t)]<Z{r,t)\Z{T,t)> . 

ti<t<tf r 

This form of path intergal representation will be the starting point of our consideration. 



5 Nonlinear representation of supergroup in Hubbard 
model 

In construction of supercoherent state we have the foUowng supermatrix which we must 
compute analyticaly: 



U — exp 



X* h+ X 

X* h- -h^ X 
\ E- X* X* -E, ) 
In this expression we have the fields of different statistics: for example, set of the fields 



(X*(^,2;,y,z),x(^,a:;,y,^)) 

are the odd grassman valued function of space-time coordinates and describe the fermionic 
degree of freedom but fields 



{E^{t, X, y, z), E+{t, X, y,z),E {t, x, y, z), h^{t, x, y, z), h'^{t, x, y,z),h {t, x, y, z)) 

describe the electro-magnetic degree of freedom, equal to two three component vectors of 
the space-time coordinates and are bosonic. In this paper we concentrate in calculating exact 
representation of exponent of the supermatrix in the coherent state. We take the dynamic 
electrical, magnetic and grassmann fields which depend on coordinates of 4-dimensional 
space-time manifold on definite the coordinates and omit {t, x, y, z) coordinates in subsiquent 
formulas. 



5 



Our general strategy will be to expand the supermatrix to N-order in fields. Them we 
can isolate and collect certain series and get some recurrent formular for general term in 
infinite series. Using this formula we can sum all terms to anytical compact representation. 
Analytical representation of the supermatrix elements will be the final point of our work. As 
a starting point we have the following exponential expression for the representation of super 
extension of the Lorentz group in the spinor representation. 

Expanding this exponent in series we can obtain first and second order in the parameters 
b and h. We see that the polynomial series on the grassmann numbers can be classified 
in grassmann order n. All the supermatrix elements can be represented as a coefficients in 
grassmann polynomials of order n, where n=0, 1,2,3,4. 



6 Matrix series for the nonlinear representation of 
supergroup 

First of all point out that the supermatrix in exponent have two submatrix: one is the odd 
grassmann matrix and the other is the even submatrix containing only the fields of type 
Ei{t, X, y, z) and hi{t, x, y, z) type. As a first step we make expansion of the exponent for 
the even matrix. Expanding in series this exponent we can obtain first, second order and 
3,4,5 order in bosonic fields Ei,i = 1,2,3. and hi,i = 1,2,3. For example, the series for 
n=0, 1,2 has the following form: 



1 
3! 






b'E, 



b^E- 




hz 
h- 





h^h- 


1 

5! 





-K 




/ 1 




Vo 

£;+ 




-E. 
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1 

+ 2! 



V 





h^h+ 


b^E, 



b'^E- 



b^E+ 





-b^Ez J 










h^h+ 





1 
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/l^ 





+ 


V 








b- ) 






(b^ 
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1 













+ 4! 























b-) 


b''E+ 


\ 







+ 






-b^E., 



J 



here we introduce the following abbriviation b = ^JE1 + E^E~ , h = \/h1 + h'^h~ 

For the odd grassmann number we have following expansion series of exponent. We write 

here two terms for the grassmann fields (x*,x) for obtaining coefficients in higher order in 

Ej and hj . 
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/ X X \ 

X* X 

X* X 

V X* X* / 



+ 



1 

2! 



/ X*^+ + XP1 X*^+ + XP3 

X*E, + x{E- + h+ + h,) 

X*E, + x{E- + h~ -h,) 

V x>2 + x^- X*P^ + XE- 





X*E+ + x{-E, + h+ + h,) 
X*E+ + xi-E, + h- -h,) 




For even order of the grassmann variables we have following series for the composite 
bosonic fields: 

f 2xx* \ 




V y 

\ 

-2E,xx* -2E,XX* 

3! -2^^XX* -2^^XX* 

-4E.XX* + {h+ + h-)xx* J 



I 4E,xx* + ih-' + h-)xx* 



+ 



+ 



1 

4! 



/ 2mmxx* AE+E^XX* \ 

-2E,{2h, + X)xx* -^E,h+xX* 

-^E,h-xx* -2E,(-2h, + X)xx* 

V -4£;-£;^XX' O 2mmxx* / 



here m — e^ + h'^,X — h'^ + h~,pm — hz + h~,pp — hz + h'^, mp — — hz, mm — h~ — hz 
and mm = m + 2E^ + E^X The expansion series for third order of the grassmann variable 
has the following form of first two terms 



/ 

E-x + {Ez + pp)x* 
E'x + {Ez + mm)x* 




m 
— r* 

4! 

^V+plX E+x*+P'iX \ 

E+x* + {~Ez+pp)x 

E+x* + {-Ez + mm)x 

E-x + p2x* E-x + p4:X* / 



\ 

here pi = Ez + pm, p2 = —Ez + pm, p3 = Ez + mp, p4 = —Ez + mp 

Expanding this series on 12 order in even and odd parameters we can obtain the analytical 
representation for the supermatrix elements. 
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7 Analytical representation of supergroup 

Having series representation for supergroup for high order in fields we can obtain the general 
analytical form for the matrix elements. Our task is to obtain exact dependence of the 
matrix element over component of fields: E^E^, E~; h^, h'^, h~ but not dependence over b 
and h. 

For example the general form of Un supermatrix elements in orders higher than 8 as a 
functions of the dynamic even and odd grassmanian fields are given by following expressions 

Uii = cosh{e) + EzSinh{e)/e+ 

(2/3 + X/2 + 2E,{2f2EE + XflEE) + 2El{2flEE + XfEE))xx* 

here the coefficients fi,i = 2,3,4 are some series in e and h variables: 
f2EE, flEE, fEE] f2h have some infinite series over e and h variables. 
Another matrix element equal to 

U21 = X*{h + {E, + /i+ + h,)h + E^ih"- + h,)f2) + xE-{h + {h"- + h,)f2) 

and we have some function for the coefficients: /4, fz, /2 
For uzi we have 

%i = x*(/4 + (i?. + h-- K)h + E,{h- - K)h) + xE-{h + + 

Last element of first column is 

W41 = E-sinh{e)/e + 2E,E-(2flEE + XfEE)xx* 

For coefficients in this expression we have: X = h'^ + h~ 
For second column we have 

ui2 = x(/4 + {E, + h- + h,)h + E,{h- + h,)f2) + X*E+{h + {h- + h,)f2) 
For 1^22 we have: 

U22 = cosh{h) + h,sinh{h)/h + E^{f2h + Xhjhh + {2h, + X)flhh){-2xx*) 
For 1*32 we have: 

= h-sinh{h)/h + E,{f2 + 2h- flhh + h-Xfhh){-2xx*) 
For M42 we have: 

U42 = X*(/4 + + h- + K)h - E,{h- + K)f2) + xE-{h + + h,)f2) 
For iii3 we have: 

1^13 = X(/4 + {E, + h+- K)U + - K)f2) + X*£;+(/3 + - /i.)/2) 
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For U23 we have: 

= h+sinh(h)/h + E,(f2 + 2h+flhh + h+Xfhh){-2xx*) 
For 1^33 we have: 

M33 = cosh{h) - h,sinh{h)/h + E,{f2h - Xhjhh + {-2h, + X)flhh){-2xx*) 
For U43 we have: 

1^43 = X*(/4 + i-E, + h+- h,)h - E,{h+ - K)f2) + xE'ih + {h-^ - h,)h) 
For we have: 

1^14 = E+sinh(e) /e + 2E,E+ (2flEE + XfEE)xx* 

For 1^24 we have: 

U24 = x(/4 + i-E, + /i+ + /i,)/3 - E,{h+ + K)f2) - X*E^{h + {h+ + /l,)/2) 
For 1^34 we have: 

K34 = x(/4 + (-^. + - h:)h - ^.(/i- - /i.)/2)+ 
X*E+{h + {h- -K)f2) 

For 1^44 we have: 

M44 = cosh{e) — EzSinh{e)/e+ 
(2/3 + Xf2 - 2E,{2f2EE + XflEE) + 2El{2flEE + XfEE))xx* 

We see later that many series in our hst arc equivalent to each other. After such selection 
between similar ones we have only some different series. 

8 Analytical representation for series 

Collecting the terms in scries expansion for a and b coefficients to 12 order we obtain for 
example the following representation for: 

1 b^ + h^ b^ + + 6^ + + h^'' + /i^ 68 + + + h%'' + 
■^^ ~ 3! ^ 5! ^ 7! ^ 9! ^ 11! ^ 

m + 
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^12 ^ ^,2^10 + ^,4^8 + /,6^6 + ;,8^4 ^ ^,10^2 ^ ^,12 

^ + 

^,14 + /,2^12 ^ ^4^10 + ^^6^8 + /^8^6 + ^^10^4 ^ ^12^2 ^ ^14 



+ 



17! 

Let us show how to sum following infinite series for b for example. We have 



52 54 56 ^8 ^10 ^2 ^4 ^6 ^8 ^10 ^12 

^+2! + 4! + 6! + 8! + 10!+ + ^^^^ + 3! + 5! + 7! + 9! + U! + 13! + ^ 



It is seen that the first series equal to 



^2 ^4 ^6 ^8 ^10 

co.M&) = l + ^ + 4y + ^ + ^ + ^ + 

and the second series equal to 

^2 ^4 ^6 ^8 ^10 ^12 

.znM&)/6 = l + - + - + - + - + - + - + 

For sum of two series we have the following expression cosh(6) + £^2sinh(6)/6 
The main series for us is the following expansion: 

1_ b^ + h^ b^ + h^b^ + b^ + + h'^b^ + b^ + + h^"" + + /i« 
~ 5! + 7! + 9! + 11! + 13! + 

bW + ^2^8 + ^4^6 + /,654 ^ ^8^2 ^ ^10 ^ 



15! 

^12 + /,2^10 ^ ^4^8 ^ ^6^6 ^ ^8^4 ^ ^10^2 ^ ^12 

17! 

^,14 + /,2^12 + /^4^10 + /.e^S + /,8^6 + ^,1054 ^ ^^12^2 + /^14 



+ 



19! 

Let us show how the summation of this series can be performed. We can make the 
summation of subpart of hole series: 

^5 ^7 ^9 ^11 ^13 J,3 

5! + 7! + 9!+Tl! + 13!+- = ""^(')-'-3! 

/i (-r + -r H r H 7 + ■■■■) = /i (smh(6) -b -) 

^7! 9! 11! 13! ^ v w 3! 5!^ 

^(^+in+l^ + --) = ^^^^^^(^)-^-^-^-^) 
Having this series representation we can rewrite expression for f 
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5! ^ 7! ^ 9! ^ n! ^ 13! 

510 + /,258 ^ ;^456 ^ /^6^4 ^ ^^852 ^ 

15! ^ 
J,12 ^ ^2J,10 ^ ^,4^8 ^ /,6^6 ^ ^8^4 ^ ^10^2 ^ ^12 



-(smh(6) - - ^) + ^/.^(smh(6) - 6 - - - -) + -h\smh{b) -b------)^ 

1^2 1^4 l3 u2 ^4 

sinh(6)/6^(l + ;^ + ;^ + •••) + - + + + -)+ 

^-5! + 6^ + F + ^ + ^-7!^ 6^^^ + 6^ + F + + - 

It is seen that the series of the type 1 + ^ + ^ + describe geometric series and 

gives the following result: 



¥ 1-^ h'-h? 

If we insert this result we obtain the representation for: 



1_ b^ + h^ 64 ^ /^2^2 ^ /^4 ^6 ^ /^2^4 ^ ^4^2 ^ ^6 58 ^ /^2^6 ^ /^4^4 ^ ^6^2 ^ ^8 

5! ^ 7! ^ 9! ^ n! ^ 13! ^ 



^10 + /,258 ^ /^4^6 ^ /^654 ^ ^^852 ^ ^10 
15! 



+ 



^ / . ^\/.5 ^'^ / b\h' b' , b\h^ b' 

— h (—b )/b h ( ) h ( ) h ... 

1,2 _ 1^2 3!^/ 1,2 _ 1^2 ^y 5!^57 52_/^2 Y\^b^b^-h^ 

sinhjb) 1 1 1,2 1,4 i 



63 62 



V 01 CI Ti ■■■v ; 



i,2_;^2 Ij2_i^2 ' y 3! 5! 71 ■■■■'i)2_i^2 



sinh(6) 1 
63 62 



+ (-sinh(/i)//i3 + 



1 



62-/^2 i,2_/^2 ■ V V // 1^2' _ 1^2 

sinh(6) sinh(/i) ^ 
~63 P~ _^ ^ 



Let us consider two series: one is / and second is /i. If we multiply / by coefficient 
and make following substitusion b— > ba, h— > ha we obtain the following series 
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+ a' ; h ; h a ; h 

5! 7! 9! 11! 

,^b^ + + + + 

« 131 + 

« ^ + 

Oj . I 



17! 

It is obvious that if we take derivative we can reduce factorial in our series for example: 

aa oa oa oa 

here we must insert b— > ab, h— > ah and a put to 1 after differentiation. 
Taking derivatives we obtain the analytical expression for : 

^ £2ff^_£i^ _ ^ giHM-^ ^ ^ cosh(?;)-c-()sh(/.) _ 

6sinh(6) — hsmh(h) 
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62 -/i2 

Series for fi have the following forms: 



■^^ ~ 4! ^ 6! ^ 8! ^ 10! ^ 12! ^ 

i,10 + /,2j,8 ^ ;,4j,6 ^ ;,6j,4 ^ ;^8^2 ^ ^10 

h 



14! 

i,12 + /,2j,10 ^ ;,4j,8 ^ ^,6^6 ^ ^,8^4 ^ ;^10j,2 ^ ^12 



16! 

fcl4 + /,2^12 ^ f^4^10 ^ ^6^8 ^ ;^8^6 ^ ^^10^4 ^ ^12^2 ^ ^14 

17! 



+ .... 



1 b^ + h-" b^ + h^b^ + h'' b^ + h%^ + h%'' + h^ b^ + h%'' + h%'' + h%^ + h^ 
•^^ ~ 3! ^ 5! ^ 7! ^ 9! ^ 11! ^ 

ftlO + /,2j,8 ^ ^,4^6 ^ ;,6^4 ^ ^8^2 ^ ^10 

h 



13! 

;,12 + /,2^10 ^ ;^4^8 ^ ;^6fc6 ^ ^^8^4 ^ ^10^2 ^ ;^12 

15! 

^,14 + /,2^12 ^ ^4^10 ^ ^6^8 ^ ^8^6 ^ ^10^4 ^ ^^12^2 ^ ^14 

17! 



+ 
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1 62 + /^2 54 ^ ^2^2 ^ ^4 56 ^ /^2^4 ^ /^4^2 ^ 58 _^ ^2^6 ^ /^4^4 ^ ^^652 ^ ^8 

■'^ ~ 2! ^ ^! ^ 6! ^ 8! ^ 10! ^ 

h 



12! 

bl2 ^ ;,2^10 ^ ;,4j,8 ^ ;,6^6 ^ ^,8^4 ^ ^10^2 ^ ^12 

W. ' 

fell + + /, 1^,10 ^ /,6,,8 ^ ;,S;,0 ^ /,10,yl ^ /,12;,2 ^ 



16! 

b^ + h^ b'^ + h^^ + h^ b'' + + + h'' b^ + h^'' + + + 
^^ = ^ + ^!~ + 4^ + 6^ + 8^ + 

+ 



10! 

fcl2 + h%10 + h^^S ^ ;^6^6 ^ ;^8^4 ^ ^10^2 ^ ^12 

12! 

i,14 + /,2^12 ^ ^4^10 ^ ^,6^8 ^ ^,8^6 ^ ^,10^4 ^ ^12^2 ^ ^14 



+ 



14! 

If we take the series for / and multiply it by 6" and take following derivative ^qq,/^ we 
obtain the series for fEE. 

1 262 + /i2 36^ + 2/i2fo2 + /i4 Ab^ + 3h%^ + 2h%'^ + h^ 

fEE = \ \ \ h 

•' 5! 7! 9! 11! 

56« + Ah%^ + Zh%'^ + 2h%'^ + /i« 

h 



13! 



15! 

For series fhh we must make multiplication of / on h? and make derivative on h? 



fhh-^ 



1 2/l2 + 62 3/^4 ^ 2fe2/i2 ^ ^4 4^6 ^ 3^,2^4 ^ 264/^2 ^ ^6 

^^^=5! + ^^ + 9^ + + 

5h^ + Ab'^h^ + Sb^h"^ + 26^/^2 + fc8 ^ 



13! 

6/iio + 562/i8 + 46^/^6 + 366/i4 + 268/i2 + ^lo 
15! 

^ ^ ^ 1 26^ + /i2 3^,4 ^ 2/i2{,2 ^ ^4 4^6 ^ 3^2^4 ^ ^h^"" + 

Jo EE = h 1 1 h 

■'^ 3! 5! 7! 9! 



13 



11! ^ 
66^° + 5hH^ + Ah%^ + Zh%'^ + 2/i«62 + /^lo 
13! 

To evaluate flEE, f2EE, flhh we multiply fEE by and make the following substi- 
tution b— > ba, h— > ha. After calculation we fix a = 1. 

It is seen that expression for flEE, f2EE, flhh are the following: 

^^^^ 1 2&2 + /i2 364 + 2/iV + /i^ 4b^ + 3h%^ + 2h%'^ + 

^^^^=4!+— 6^ + 8^ + m + 

56^ + Ah^^ + 3h%'^ + 2h%'^ + ^ 



12! 



14! 

d'^fEE 



+ 



Repeating all operation we can obtain for f2EE — (— ^^)a=i 

, ^ ^ 1 262 + + 2//-/r + /i^ 4// + 3//'-6i + 2/; ^/r + /^.^ 

fgEE = \ \ \ h 

■'^ 3! 5! 7! 9! 

56« + Ah%^ + 3/1^6^ + 2/1^62 + 

h 



11! 

66^° + hh%^ + 4/^^66 ^ 3/^6^4 ^ 2/i862 ^ ^lo 



13! 

and for flhh = {^)a=i and for f2hh = (^)„=i 

1 2/1^ + 62 3/^4 ^ 2fe2/i2 ^ ^4 4^6 ^ 3^2^4 ^ 264/^2 ^ ^6 

^^^^=4!+^^ + 8^ + 10^ + 

5/i« + 462/i6 + 3^4/^4 + 26^/^2 + 6^ ^ 



12! 

6/i^° + + 46^/1^ + 366/i^ + 26^/12 ^ ^lo 



14! 



+ 



1 3/^2 + 62 5/^4 ^ 3^2/^2 ^ J,4 7/^6 ^ 5J,2;^4 ^ 3J,4^2 ^ ^6 

^'^=3!+^^ + v. + 9^ + 

9h^ + n^h^ + 56^/1^ + 36^/^2 + ^,8 

h 



11! 

ll/i^" + 962/^8 + 76^/1^ + hb'^h* + 36^/12 + 



+ 



13! 

In analytical form we have for f2h — 2f2hh — f2 

Series for DEDEf equal to derivative of fEE on 6^. It is seen if we compare series for 
fEE and DEDEf: DEDEf = 

14 



^^^^^ 2 + 2/1^ 12&4 + 6/i2&2 + 206^ + 12^ + 6/^4^,2^2/16 
DEDEf=- + ^^ + + + 

306« + 20hH^ + I2h%^ + 6fe%2 + 2h^ ^ 



15! 

We can obtain the following representation for DEDEfl : DEDEfl = 

2 662 + 2/i2 126^ + 6/1252 + 2/i4 206^ + 12/126^ + 6/1^6^ + 2/1^ 

DEDEfl = — H ; ^ ; \ ; h 

■' 6! 8! 10! 12! 

306« + 20/12^6 + 12/i4fo4 + 6/^6^,2 + 2/i8 

w. + 

Comparing series for DEDhf2 and series for /2 we can obtain DEDhf2 — ^^^f^^^^a-^ 

^ ^ ^ , 1 262 + 2/l2 3^4 + 4^2^2 + 3/^4 4^6 + g/^2 + g/^4j,2 + 4/^6 

i^^W2 = ^ + ^^ + + 13^ + 

56^ + 8/i266 + 9/^4^4 + 8/^6^2 + 5/^8 

h 



15! 

66^° + 10h%^ + I2h%^ + I2h%^ + 10/1^62 + e/i^o 



+ 



17! 

and for DEDhfS we can obtain the following representation: DEDhfS — d{tfi)dlh^) 



^^^^^3 = ^ + ^^+ m + m + 

56^ + 8//-//' + 9]iH>^ + 8//*Vr + r)/,8 

+ 



14! 

66^° + lOhH^ + I2h%^ + I2h%'^ + 10/i«62 + 6/1^° 



+ 



16! 

Summing our calculation let's write all analytical formula for the function in the matrix 
elements: 

2 sinh(/i)63 + /i (62 _ /i2) cosh(6)6 + {h^ - Zh'^h) sinh(6) 
~ 263/1 (62 -/i2) 2 ' 

-26cosh(6) + 26cosh(/i) + (62 - /i2) sinh(6) _ 
~ 26 (62 -/i2)2 ' 

h{h'^- h?) cosh(6) - (62 + /i2) sinh(6) + 26/isinh(/i) 
" 26 (62 -/i2)2 

/isinh(6) — 6sinh(/i) cosh(6) — cosh(/i) 

" 63/i - 6/i3 ' " 62^r7^2 ; 
6sinh(6)-/.sinh(/.) ^Hffi) _ Elffi) 

62 -/i2 ' 62 -/i2 

15 



9 Conclusion 



We have calculated the exact representation of the supergroup as well as the supercoherent 
state in the Hubbard model. These constructions naturally appear in the strongly correlated 
electronic systems in the case of introducing atomic bases in limit of large on-site Hubbard 
repulsion. The dynamical supergroup which operates in a local superbundle determined by 
any on-site cigcnfunction gives us the wave function in the form of a superspinor. This super- 
spinor describes a local supercoordinates frame in the curved supermanifold. The operator 
spinor part acting in tangent and cotangent bandies of this supermanifold in supergroup 
can be reformulated in the terms of the atomic Hubbard operators. Next step lies in the 
calculation of effective functional of Hubbard model. 
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